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Abstract 

The rise of quantum information science has opened up a new venue for applica¬ 
tions of the geometric phase (GP), as well as triggered new insights into its physical, 
mathematical, and conceptual nature. Here, we review this development by focusing 
on three main themes: the use of GPs to perform robust quantum computation, the 
development of GP concepts for mixed quantum states, and the discovery of a new 
type of topological phases for entangled quantum systems. We delineate the theoretical 
development as well as describe recent experiments related to GPs in the context of 
quantum information. 
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INTRODUCTION 


A cat held upside down can fall on its feet although its angular momentum vanished. An 
amoeba can swim in a completely reversibel environment®. How do they manage? 

The mechanism used by these cats and amoebae is the same as that used by a quantum 
state that returns to itself, but picks up a phase factor of geometric origin. The mechanism is 
called holonomy or geometric phase (GP) and describes the twisting of a physical quantity, 
such as the orientation of the cat or the phase of the quantum state, due to the curved 
geometry of the physical space where the change takes place®. This space is the space of 
shapes in the case of cats, while it is the quantum state space or projective Hilbert space in 
the case quantum states. 

If the evolution of a quantum system is induced by slowly changing some external pa¬ 
rameters around a loop, the state approximately performs a loop giving rise to an adiabatic 
GP. This kind of GP was first discovered in the field of quantum chemistry by Longuet- 
Higgins and coworkers®® and Stone®. It is the underlying mechanism behind the molecular 
Aharonov-Bohm effect®, which takes the form of a topological phase shift when a molecule 
reshapes slowly around a conical intersection point. 

Later, Berry® demonstrated that the GP is a general consequence of the geometrical 
structure of the space of parameters that drive cyclic adiabatic evolution of energetically 
non-degenerate quantum states; a scenario that can be realized in numerous systems, such 
as photons®, NMR®, neutrons®, Jahn-Teller molecules®, condensed matter system^®, and 
cold atoms®. Simon® demonstrated that Berry’s GP can be understood as the holonomy 
in a Hermitian line bundle. This result has established the geometrical nature of the phase. 

GPs have been generalized beyond Berry’s framework. Wilczek and Zee®, still within 
the context of adiabatic evolution, pointed out that energetically degenerate states may 
aquire matrix-valued GPs when slowly changing parameters trace out a loop in parameter 
space. The matrix nature of these phases make them potentially non-commuting; a fact 
that has motivated the term non-Abelian GP or quantum holonomy for this kind of phase 
effects. Aharonov and Anandan® removed the restriction of adiabatic evolution and intro¬ 
duced the concept of non-adiabatic GP in arbitrary cyclic quantum evolution. Subsequently, 
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Anandanl22l proposed non-adiabatic non-Abelian GPs of general quantum systems. Based 
on Panchartnam’s work*^ on polarization of classical light fields, Samuel and Bhandari^^ 
defined GPs for open paths in the underlying state space. 

The rise of quantum information science has opened up a new venue for application of the 
GP, namely, to use it as a tool for robust quantum information processing^^SM] devel¬ 

opment has been triggered by the work of Zanardi and coworkers^^^® on holonomic quantum 
computation, i.e., the idea to use non-Abelian Wilczek-Zee GPs to implement quantum gates; 
the logical transformations that build up a circuit-based quantum computation. Such im¬ 
plementations are believed to be useful to reach the error threshold, below which quantum 
computation with faulty gates can be performed. The basic reasoning behind the conjectured 
robustness is that GP is a global feature of quantum evolution being resilient to errors, such 
as parameter noise and environment-induced decoherence, which are picked up locally along 
the path in state space. 

Gonversely, tools and concepts developed in quantum information theory have been used 
to broaden the concept of GP itself. Based on early work by Uhlmann^ on holonomy along 
paths of density operators, new concepts of GPs for statistical mixtures of wave functions 
have been developed^^®^. In relation to these mixed state geometric phases, GPs of quantum 
systems that interact with a quantum-mechanical environment have been examined in the 
contexts of quantum jump^^^*^, quantum maps^^^^, stochastic unravellings^^^^, and the 
adiabatic approximationSS®. Studies of GPs of two or more quantum degrees of freedom 
have led to the discovery of entanglement-induced topological phase^^. This new type of 
phases have been shown to be useful to characterize quantum entanglement. 

The objective of this review is to describe the recent merging of ideas in quantum in¬ 
formation science and the held of GP. The basic theory of Abelian and non-Abelian GPs 
are outlined in the next section. Thereafter follows the core of this review, namely, a de¬ 
scription and overview of various applications of GPs in quantum information. This part is 
focused on three main themes: quantum computation, mixed quantum states, and quantum 
entanglement. The paper ends with a concluding section, which in particular contains some 
pertinent issues to examine in the future in this research area. 


3 






DIFFERENT FLAVORS OF GEOMETRIC PHASE 


The general structure of the quantum geometric phase (GP) is the removal of accumulated 
local phase changes from the global phase acquired in some evolution of a quantum system®! 

GP = Global phase — Local phase changes. (1) 

The resulting GP factor can be an Abelian phase factor, denoted as or a non-Abelian 

unitary matrix, denoted as G[C], depending on the context. The path C may reside in 
different kinds of spaces, such as projective Hilbert space in the case of pure states, the 
space of subspaces in the case of quantum gates, and the space of density operators in the 
case of mixtures of pure states. In this section, we review the basic theory of these Abelian 
and non-Abelian GPs by illustrating them in the physical context of pure state evolution. 

Abelian GPs 

A well-known feature of quantum evolution driven by a time-independent Hamiltonian H 
is that a stationary pure quantum state, as represented by a vector |?t,) in Hilbert space of 
the system, picks up a phase factor determined by the average energy {n\H\n) = E and the 
elapsed time f, i.e., 

IV’(O)) = \n) iV’(t)) = (2) 

This simple fact is not true when the state itself evolves in time as the phase then picks up 
information of the geometry of the path in state space in addition to the system’s energy. 
The origin of this extra phase contribution, being the pure state GP, is the curvature of state 
space. 

To understand the appearance of a GP shift accompanying state changes, consider a 
superposition of two stationary states \n) and \m) of the form 

\^jJ{0)) = a\n) + b\m). (3) 

Let En and Em be the corresponding energies assumed to be non-degenerate and ordered as 
Em > En- The coefficients a and b are non-zero complex numbers whose values determine 
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the initial state of the system, such that |ap and | 6 p are the probabilities to hnd the energies 
En and respectively. Linearity of the Schrodinger equation implies that |t/’(0)) evolves 
into 




(4) 


This describes a non-trivial cyclic evolution of the quantum state as can be seen by evaluating 
\'ip{t)) at time r = 2nh/{Em — En): 

|^/’(r)) = {a\n) 


r is thus the period of the cyclic evolution that corresponds to the loop C : [0,t] 3 t ^ 
|'^(t))('^(t)| in state space. 

The resulting overall phase / = —27rEn/(Em — En) can be divided into a sum of two parts. 
The hrst part is the phase induced by the average energy {'ip(t)\E[\'ip(t)) = \a\^En + \b\‘^Em 
and reads 


O-TT 

■5=(lorsn+itrs™) (6) 

in analogy with the above phase shift Et/h for a stationary state. <5 is clearly energy- 
dependent and therefore a dynamical phase associated with the evolution. The remainder 
7 = / — 5 reads 


7 = 


27rE„ 


+ 


2tt 


{\a\^En + \b\^E^)=27r\b\^, 


(7) 


where we have used normalization |ap -|- | 6 p = 1. Remarkably, 7 is independent of the 
energies of the system; in fact, 7 is the GP <h[C] associated with the loop C. 

If we remove the second phase factor instead of the phase factor then 

the overall phase / = —2'KEm/{Em — En) and the GP 7 = —27r|ap 7 ^ 7 . By using the 
normalization of |'^(0)), we hnd 7 — 7 = 27i, which implies = e^'^. In other words, the GP 
factor is the same for the two arbitary phase choices. Indeed, it is not difficult to prove that 
it is invariant under any phase choice of the evolving state, which makes the GP of a pure 
state an invariant under the Abelian group U(l), i.e., the group of phase transformations. 
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The independence of energy and local phase choice suggest that the Abelian GP can be 
dehned without reference to the underlying dynamical mechanism that drives the evolution. 
This can be seen by using the concept of relative phasd^ between state vectors. Let \i/j{ti)) 
and \‘ip{t 2 )) be two non-orthogonal vectors along a continuous curve G : t G [0, r] — )■ |'0(t)) in 
Hilbert space. The phase of \ijj{t 2 )) relative |'0(G)) is simply the phase of the scalar product 
(' 0 (ti)|' 0 (t 2 ))- The relative phase concept can be used to dehne 

Global phase = arg('0(O)|'0(r)), ( 8 ) 


being well-dehned if ('0(O)|'0(r)) 7 ^ 0 (non-orthognal end-points), and 

nr pT 

Local phase changes = lim / a.Tg{'ijj{t)\‘ip{t + 6t)) =—i / {'ip{t)\‘ip{t))dt, (9) 

^ < 5*^0 Jo Jo 

where we have used Taylor expansion -|- 6t)) 1 -f {-ifj{t)\'ip{t))6t ^ and 

the fact that (' 0 (f)|'^(t)) is purely imaginary due to normalization; 1 = (' 0 (t)|' 0 (f)) ^ 0 = 
{'ip{t)\'ijj{t)) + {il>{t)\‘ip{t)) = 2Re('0(t)|'0(t)). By combining Eqs. (|^ and (|^ with Eq. Q, we 

obtain®! 

d>[C] = arg(?/>(0)|V’(r))+i / {^p{t)\ 1 p{t))dt. (10) 

This GP is invariant under local phase changes |t/’(t)) —)■ e*"*'*^|'^(t)) and invariant under 
reparametrizations t —)■ r(t) such that r > 0,Vt G [0,r]. Thus, GP is a property of the path 
C in state space, dehned by the projection \'ip(t)) —)■ \'ip(t)){'ip(t)\ of the path C in Hilbert 
space. Note that this dehnition of GP makes no reference to cyclic evolution and therefore 
applies to any evolution of a pure quantum state connecting two non-orthogonal end-points 
|?/^( 0 )) and 

There are three ‘canonical’ choices of local phase: 


(i) If I'lpit)) is a solution of the Schrodinger equation, i.e., satishes ihl'ipit)) 
then we obtain from Eq. (10) the expression®! 




d>[C] = arg(V’( 0 )|^/>(r)) - \ ( 11 ) 

where H(t) is the Hamiltonian of the system. Thus, the GP can be understood as the 
removal of the accumulated phase induced by the average energy {'ip(t)\H(t)\'ip(t)), i.e.. 
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the dynamical phase shift of the system, from the acquired global phase, just as in the 
example of the two superposed stationary states discussed above. 


(ii) If 


(V’(t) \ij{t + St)) > 0 ^ \'ip{t)) = 0, te [0, r) 


no phase changes are acquired along the path, and we obtain 


( 12 ) 


t[C] = arg(i/i(0)|i/i(T)). 


(13) 


When Eq. (12) is satished, \iS{t)) is said to be parallel transported. In this case, the 
GP coincides with the global phase associated with the evolution. Parallel transport is 
particularly useful in experiments as it makes the GP directly accessible by observing 
the global phase^^^. 


(hi) The phase choice |A(f)) = |'0(t)) is a gauge invariant reference section^S^I^ 

as |A(f)) is unchanged under local phase transformations of \'ijj{t)). (To be precise, 
|A(f)) —)■ e*"^°)|A(f)) when \7p) —)■ We note that arg(A(0)|A(r)) vanishes so 

that 


d>[C]=* [\x{t)\Xit))dt, (14) 

Jo 

which is the familiar form of the geometric phase being the closed path integral of the 
effective ‘vector potential’ i{X\dX). This perspective was adopted in Berry’s original 
work on GP®. 


Note that although the three canonical phase choices focus on different aspects of quantum 
evolution, they differ only by local phase choices and therefore all give the same numerical 
result for the GP. 


Example: GP of a qubit 

A qubit is a quantum system that can be described by two states |0) and |1). Physical 
examples of qubits are abundant: it could be the polarization of a photon, the spin of an 
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electron, neutron, or proton, two isolated states of an atom or ion, two charge state of a 
superconducting island, and so on. An arbitrary qubit state \'ip) can be written as 

IV’) = cos ^|0) + e*'^sin ^|1). (15) 

The two parameters 9 and 0 may be interpreted as polar angles of a unit sphere, called the 
Bloch sphere, with |0) and |1) projecting onto the north {9 = 0) and south poles {9 = |), 
respectively, see Figure 
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Figure 1; Bloch sphere representing pure qubit states with the projection of the basis states 
|0) and |1) at the poles. The angles 9 and 0 determine the relative weight and phase shift, 
respectively, between |0) and |1). A loop C enclosing the solid angle on the Bloch sphere 
is associated with the GP — 


Now, by assuming that the qubit traces out a loop C : f G [0,r] —>■ ( 6 ^ 4 , 0i) on the Bloch 




sphere, i.e., Or = Oq and (j)r = 0o, we obtain 

<h[C] = — '/(I — cos 9)d(() = — Q, (16) 

2 ic 2 

f2 being the solid angle enclosed on the Bloch sphere by C. By nsing Stokes theorem, we 
obtain 

<1)[C] = —^ <£ smOdOdcj) (17) 

2 Js 

where S is the snrface on the Bloch sphere enclosed by C. This is eqnivalent to the flnx 
throngh 5 of a magnetic monopole of strength — | sitting at the origin. 

Non-Abelian GPs 

Following Berry’s work®' on Abelian GP factors accomanying adiabatic changes, Wilczek 
and Zee^ pointed ont that non-Abelian gange strnctnres appear when more than one state 
is considered. This tnrns ont to be relevant for instance in adiabatic evolntion where several 
energy eigenstates are degenerate over some region of parameter space. 

The key ingredient of the non-Abelian generalization is the evolntion of a snbspace of 
the fnll state space. Let t G [0, r] —)■ St trace ont a loop C (i.e., Sr = Sq) of iP-dimensional 
snbspaces of an A^-dimensional Hilbert space. Let {\'ipk{t))}k=i span St- Since C is a loop, the 
overlap matrix Ufci(0,r) = ('0fc(O)|'0«('^)) is nnitary and corresponds to the overall ‘phase’ of 
the evolntion. By applying the same reasoning for two inhnitesimally close instances t and 
t+6t, we obtain the local ‘phase shift’ to hrst order in 6t. Here, Afcz(t) = i{^jJk{t)\^jJl{t)). 

The accnmnlated local phase shift along the loop is Te^-io where T stands for time 

ordering. Thns, we obtain the non-Abelian GP; 

V[C] = lJrS^oMt)dt^ (18) 

where U = U(0, r). Note that for the special case where K = 1, U[C] rednces to the complex 
nnmber 

which we recognize as the Abelian GP factor discnssed above. 
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The subspace need not trace out a loop to acquire a non-Abelian GP. For such an open 
path C, the overlap matrix Mfc;(0,r) = is no longer unitary, but admits a 

unique polar decomposition M = |M|U provided |M| > 0. We therefore have 

= |M|Ure*^o"A(t)rft = |m|U[C], (20) 

which dehnes the non-Abelian GP® 

U[C] = \M\-^Mre^^o Mt)dt ('21) 

valid for any evolution connecting two fully overlapping A'-dimensional subspaces Sq and 
Sr of the A^-dimensional Hilbert space. Note that the existence of the inverse |M|“^ is a 
consequence of the condition |M| > 0 that is required for a unique polar decomposition of 
the overlap matrix. 

The matrix-valued GPs may as well be viewed as operators: 

K 

U[C] = Y. Ut,lC||Vi(r)>(V,(0)|, (22) 

k,l=l 

with initial subspace as domain and hnal subspace as image (these subspaces coincide if C 
is a loop). The matrix and operator forms of the non-Abelian GPs are fully equivalent. 

The non-Abelian GP transforms as U[C] —)■ V(0)U[C]Vl(0) under a unitary change of 
basis |t/’A:(t)) —)■ Yhi along the moving subspace St. This implies that the Wilson 

loop TrU[C] is unchanged under local basis changes. Thus, the non-Abelian GP transforms 
gauge covariantly and is therefore a property of the path C in the space of subspaces®. 

Example: Geometric transformation of a qubit 

A technical complication when calculating the non-Abelian GP is the need to evaluate a 
time-ordered product. This may be dealt with by performing a pair of loops resulting in 
non-commuting GPs, therefore demonstrating the desired non-Abelian feature, but along 
each of these loops no time-ordering is needed. Here, we illustrate this idea in the case of 
a two-dimensional subspace 5, encoding a single qubit, that evolves in a three-dimensional 
Hilbert space "H. 
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Let T-L = Span{|0), |1), |a)}. Consider paths of two-dimensional snbspaces C ^ St = 
Span{|'^(t)), snch that 5o is the qnbit snbspace. Thns, Sq = Span{|0), |1)}. 

Let ns hrst consider the choice 

l^iW) = |o), 

= cosyll)+e*’'*sin^|a). (23) 

A loop Cl that starts and ends at 5o is implemented by choosing = kq = 0 and ? 7 r = ho- 
The GP associated with Ci is 

U[Ci] = (^24) 

The phase shift can be interpreted as the solid angle enclosed on the Bloch sphere associated 
with the snbspace Span{|l), |a)} and parametrized by the spherical polar angles k,t]. 

Now, consider instead the choice 

|'02(i)) = -sin?7t|0)cos?7t|l), 

= cosKt cosht|0)-|-cosKtsinr/fll) — sinKt|a), (25) 

which span a snbspace that starts and ends at Sq provided k,- = kq = 0 and hr = ho- The 
GP associated with this loop C 2 is 

(26) 

where ay = —i|0)(l| -|-i|l)(0|. Jnst as G[Ci], the nnitary qnbit transformation U[C 2 ] has a 
geometric interpretation in terms of a solid angle bnt with a different meaning: it is the solid 
angle enclosed on a parameter sphere that dehnes all possible choices of two-dimensional 
snbspaces with real-valned coefficients with respect to the given reference basis {| 0 ), | 1 ), |a)} 
of the fnll Hilbert space. 

For non-trivial loops Ci and C 2 the resnlting U[C]\ and t/[C 2 ] sue non-commnting thereby 
demonstrating the desired non-Abelian featnre of these GPs. The above qnbit GPs were 
discovered theoretically by Unanyan et alJ^ in the context of atoms controlled by stimnlated 
Raman adiabatic passage (STIRAP) techniqnes. 
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QUANTUM COMPUTATION 


A classical computer works by manipulating a set of bits according to some algorithm to 
produce a desired result. Analogously, a quantum computer manipulates qubits by applying 
a certain combination of logical operations, quantum gates, forming a quantum circuit. 
GPs have been proposed to be useful in realizing quantum gates. One of the main reasons 
for considering GP-based quantum gates is that they are believed to be more robust than 
traditional dynamical gates, and therefore useful in order to reach the error threshold^, 
below which quantum computation with faulty gates can be performed by means of error 
correction protocol^^^®. 

A key ingredient of circuit-based quantum computation is the concept of universality, 
which is the ability to perform an arbitrary unitary transformation on a certain set of qubits. 
It has been shown that universality can be achieved by applying arbitrary single-qubit op¬ 
erations and at least one entangling two-qubit gate^. An arbitrary single-qubit operation 
requires the ability to perform non-commuting gates. Geometric quantum computation 
(GQG) is the idea to use GPs to implement such a universal set of one- and two-qubit gates. 

There is basically two different methods to achieve universality by using GPs. Either, 
one realizes geometric phase shift gates with respect to different based^®®, or one realizes 
non-commuting gates by using non-Abelian GPs for K > 1 dimensional subspace^^^^. Each 
of these methods can be implemented by employing adiabatic or non-adiabatic evolution. 

GQC based on Abelian GPs 

Geometric phase shift gates take the form \k) —)■ with fk being the GP of the com¬ 

putational state \k). At hrst sight, such gates are not sufficient for universality since they 
all commute. To resolve this, one may instead implement GPs with respect to different 
bases®®. The resulting phase shift gates become non-commuting and therefore potentially 
universal. 

To demonstrate this idea, consider a single-qubit \k), k = 0,1. Assume the orthogonal 
states |-1-) = cos ||0) -|-e*'^sin ||1) and |—) = —e“*'^sin ||0) -|-cos ||1) aquire purely GP shifts 
=pr2/2 after cyclic evolution (G being the solid angle enclosed by the Bloch vector). Thus, 
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|±) —>■ defining the geometric gate transformation^SI] 

U{n,n) = e-^^^^\+){+\ + , (27) 

where cr = {(Jx,(Jy,az) are the Panli operators and n = (sin 6* cos 0, sin 6* sin 0, cos 0). Pro¬ 
vided the solid angle and n can be varied independently, U{Q,n) is an arbitrary SU(2) 
transformation, i.e., a universal one-qubit gate. 

The dynamical phase contributions of the |±) states can be eliminated either by employ¬ 
ing refocusing technique^^^®, rotating driving fields with fine-tuned parameterd ^^ * ^^ * ^^ , or by 
driving the qubit along geodesics on the Bloch sphere by using composite pulsed^. These 
techniques result in non-commuting gates solely dependent on the GPs of the cyclic states. 

Physical implementation: NMR qubits 

Quantum computation based on nuclear magnetic resonance (NMR) technique manipulates 
qubits encoded in nuclear spin-^. Such a qubit evolves due to the Zeeman interaction pi ■ B, 
fi being the magnetic moment of the nuclei, I the spin of the nuclei, and B the applied 
magnetic field. 

The nuclear spin states | t; n) cind | 4-;n), n being the quantization axis, acquire GPs 
—Vt/2 and -l-G/2, respectively, by adiabatically turning the magnetic field around a loop 
C that starts and ends along n and encloses the solid angle G. The dynamical phases are 
eliminated by refocusing, whereby the spin is taken around C twice, the second time exactly 
retracing the first evolution but in opposite dirrection, and performing fast spin flips (tt 
transformations) immediately before and after the second loop. Thus, the spin undergoes 
the evolution C —)■ vr —)■ C~^ —>■ tt. While the GPs of the qubit spin states | t; n) and | 0; n) 
add up after performing the spin echo sequence, the dynamical phases exactly cancel leaving 
a purely geometric phase shift gate 

G(2G,n) = (28) 

which is a universal one-qubit operation as noted above. A key point of the refocusing 
technique is that it makes the scheme insensitive to spatial variations of the magnetic held 
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strength so that all spin qubits in the sample undergo the same purely geometric transfor¬ 
mation, determined by the same solid angle swept by the direction of the slowly varying 
magnetic field^. 

Next, a conditional phase shift gate acting on two nuclear spin qubits can be realized 
by utilizing the standard NMR uniaxial spin-spin interaction Jlz^a ® Iz-b in addition to the 
Zeeman term {fia^a + hfelfe) • B of the two spins a and b. The interaction term effectively shifts 
the ^-component of the magnetic held by J/ + Pb) or — J/ (pa -|- fib) depending on whether 
the spins are parallel or anti-parallel. By choosing the initial quantization axis of the two 
spins along the ^-direction, the sequence C —)■ tTq —)■ C~^ —)■ tt;, —)■ C —)• vr^ —)■ C~^ —)■ Tib, 

where tt^ and Tib are spin hips applied selectively to spins a and b, respectively, and C is an 
adiabatic loop of the magnetic held, results in the conditional phase shift gate^ 

t/(A7) = e^^^^{\taUz){taUz\ + \Ub,z){Uib,z\) 

+e taift; 2;)(taib; 2;| + I iatb) ^)(4'atb; ^1)- (29) 

Here, 

Ay = 7+-7-, (30) 

where 7 ± = Vt±/2 with being the solid angles enclosed by the ehective magnetic helds 
{Bx, By, Bz ± J/ {fia + /ife)) • 

Note that U (Ay) is entangling as it cannot be written as a product of unitary operators 
acting locally on each nuclear spin. Thus, f/(2ff,n) and [/(Ay) form a universal set of 
one- and two-qubit gates that can be used to perform any quantum computation by purely 
geometric means. 

Experiments 

The conditional geometric phase shift gate U (Ay) has been demonstrated in NMR by Jones 
et al^. The two qubits were taken as the weakly interacting nuclear spins. The desired 
geometric gate was realized by adiabatically sweeping the magnetic held around a loop C 
starting and ending along the z axis. 

Another class of experiments related to GQC is based on the idea that the dynamical 
phase can in some cases be proportional to the geometric phasd^. This avoids the need 
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to eliminate the dynamical phase. Since such gates have non-zero dynamical phase, they 
are not geometric gates in a strict sense; nevertheless they share the potential robustness of 
standard geometric gates. The term ‘unconventional GQC’^ has been coined for this type 
of gates. 

InlSSl, a conditional vr-phase two-qubit unconventional geometric gate combined with 
single-qubit rotations to produce a Bell state with 97% fidelity of two trapped ions has 
been demonstrated. A universal set of one- and two-qubit unconventional geometric gates 
with average fidelities 97-99% (one-qubit gates) and 93% (two-qubit gate) has been demon¬ 
strated in NMR^. These high fidelities indicate the potential usefulness of Abelian geometric 
phases for robust quantum computation. 

Holonomic quantum computation: GQC based on non-Abelian GPs 

Non-Abelian GPs are matrix-valued and can therefore be non-commuting. By finding meth¬ 
ods to eliminate the dynamical phases, all-geometric universal quantum computation can 
thus be implemented by using non-Abelian GPs. 

In the standard scheme of non-Abelian GQCPHSZl^ 

the dynamical phase is eliminated 
by utilizing adiabatic evolution of energetically degenerate subspaces. In this way, matrix¬ 
valued Wilczek-Zee phases^ are implemented that can be used to realize a universal set of 
quantum gates^^®^. A difficulty in realizing these gates is the long run time associated with 
adiabatic evolution. In other words, adiabatic non-Abelian geometric gates operate slowly 
compared to the dynamical time scale, which make them vulnerable to open system effects 
and parameter fluctuations that may lead to loss of coherence. 

To overcome the problem with the long run-time, GQG based on non-adiabatic non- 
Abelian GPs has been proposed^ and further developed^^®^. The non-adiabatic scheme 
can be performed at high speed and involves less parameters to control experimentally than 
in the standard adiabatic scheme. 

In the following, we describe the basic idea of adiabatic and non-adiabatic non-Abelian 
GQG. We describe the schemes in terms of trapped atoms or ions, which provide a versatile 
tool to perform quantum information processing by geometric means. Qubits encoded in the 
energy levels of the atoms or ions can be manipulated by external fields and the trapping 
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potential to realize universal non-adiabatic non-Abelian GPs. 


Non-Abelian GQC: Adiabatic case 


Adiabatic non-Abelian GQC can be implemented in an atomic or ionic system by controlling 
transitions between four energy levels by three independent laser pulses forming a tripod 
configuration. This scheme, hrst proposed for trapped ions or atoms^^^^ (see alscP^^, and 
later for superconducting qubit^^ and quantum dots^, has become the standard one to 
perform GQC by means of adiabatic evolution. Other schemes not based on the tripod 
conhguration can be found 

The tripod shown in Figure]^ consists of four atomic or ionic energy levels |0), |1), |a) 
and |e) of the ‘bare’ Hamiltonian Hq = — h(/eo|0)(0| -|-/ei|l)(l| + fea\o){a\) controlled by 
suitable external helds (by putting the ‘bare’ energy of the excited state to zero). The states 
| 0 ), | 1 ), |a) couple to the ‘excited’ state |e) by applying three oscillating electric held pulses 
Ej(t) = ejgjit/T) cos{fjt), j = 0,1, a, where the g/s describe slowly varying shapes and 
relative phases of the pulses, and the e^-’s describe the polarization of the laser beams. By 
tuning the oscillation frequencies fj so that detunings satisfy fj — f^j = A, the Hamiltonian 
in the interaction picture reads 

= -hA\e){e\ + h{u:o{t/T)\e){Q\+u:i{t/T)\e){l\+u:a{t/T)\e){a\) + h..c., (31) 


where ujj = {e\ii ■ ej\j)gj/{ 2 h)^ /j, being the electric dipole moment operator, and we have ne¬ 
glected rapidly oscillating counter-rotating terms (rotating wave approximation 

(RWA)). 

The interaction Hamiltonian has two degenerate dark energy eigenstates: 

\Di;9,(l)) = — sin(pe**''^‘^“‘^^^|0)-f cos(pe*^‘^‘^“‘^^^|l), 

|-D 2 ; 6 ', 0 ) = cos'd (cos<pe**''^‘^“'^^^|0) -I- sin<pe**^'^‘^“'^^^|l)) — sin'd|a), (32) 


where uq = cu sin'd cos , ui = u sin'd sin , and Ua = ui cos 'de*'^^. These dark eigen¬ 
states coincide with the realization of the non-commuting one-qubit operations U[Ci] and 


U[C 2 ] in Eqs. (24) and (26) by identifying the adiabatic parameters (■d, (p, Rs — Si, S 3 — S 2 ) 
with {k/2,7i/ 2,0, —rj) and (k, 77 , 0 , 6 ), respectively. These two gates are realized by varying 
the parameters k^tj in the T —)■ 00 adiabatic limit. 
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Figure 2: Tripod setting consisting of three energy levels |0), |1), and |a) coupled to an 
excited state |e) at zero energy by laser pulses with the same detuning. Adiabatic non- 
Abelian geometric gates can be implemented by slowly varying the complex-valued laser 
parameters ujj, j = 0,1, a, around a loop in parameter space. The parameters start and end 
so that the two-dimensional dark subspace coincide with the qubit subspace Span{|0), |1)}. 

The above U[Ci] and U[C 2 ] form a universal set together with any entangling geometric 
two-qubit gate. For trapped ions, such a gate can be performed by utilizing the Sprensen- 
Mplmer setting^, resulting in the Hamiltonian^ 

Hf\t) = a; sin —e*'’|ee)(ll|-|-cos —|ee)(aa|-|-h.c. j , (33) 

which has a single dark state \D] k, rj) = cos f |11) + e'^^sin ||aa) that picks up the GP 

IJ[C-^ — g-Tl>(ll|/{l-cosK)dr; 

by adiabatically changing k and r] around a loop in prameter space. U[Ci], U[C-^, and U[C^ 
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form a universal set of geometric quantum gates that can be used to build any quantum 
computation by purely geometric means. 

Non-Abelian GQC: Non-adiabatic case 

Non-adiabatic non-Abelian GQC^^ in an atomic or ionic system requires control of three 
energy levels by suitable external laser helds forming a A-system. By following the same 
procedure as above, we obtain the Hamiltonian 

= a;o(t)|e)(0|-l-a;i(t)|e)(l|-Fh.c. (35) 

by employing the RWA in the interaction picture and by requiring zero detuning. A key 
difference compared to the adiabatic scheme is that the complex-valued Rabi frequencies 
Uo{t),Ui{t) are allowed to be varied at any speed. 

A non-Abelian geometric one-qubit gate U[C] acting on the qubit subspace is imple¬ 
mented by choosing laser held pulses of duration r such that cooit)/coiit) = — tan(6'/2)e®'^ 
is time independent and satisfy the criterion + \uji{t)\^dt = tt. Here, the latter 

condition assures that the qubit subspace Span{|0), |1)} undergoes a cyclic evolution dehning 
a loop C in the space of two-dimensional subspaces of the three-dimensional Hilbert space 
Span{|0), |1), |e)}; the former guarantees that the dynamical phases vanish for the full pulse 
duration, which implies that the gate is fully determined by C. Explicitly, 

U[C] = n ■ cr, (36) 

where cr are the Pauli operators acting on the qubit subspace. This unitary corresponds to 
a 180° rotation of the qubit around the ‘direction’ n = (sin 0 cos (/>, sin 0 sin 0, cos 6^) dehned 
by relative phase and amplitude of the pair of laser helds. Performing sequentially two such 
gates with different laser parameters dehning directions n and n' yields 

U(C')U[C] = n' ■ n -|- icT ■ (n' X n), (37) 

which is an arbitrary SU(2) operation, i.e., a universal one-qubit gate. Geometrically, the 
gate can be visualized as a rotation of the qubit by an angle —2 arccos (n' ■ n) around the 
direction n' x n. 
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Figure 3: A setting consisting of two energy levels |0) and |1) coupled to an excited state |e) 
at zero energy by a pair of zero-detuned laser pulses. Non-adiabatic non-Abelian geometric 
gates acting on a qubit encoded in the subspace Span{|0), |1)} can be implemented by varying 
the complex-valued laser parameters cuj, j = 0,1,a, but keeping ojo(t)/oji(t) constant over 
each pulse pair chosen to satisfy the tt pulse criterion = vr. 


Similar to the adiabatic realization, the universal set is completed by adding a geometric 
two-qubit gate in the Sprensen-Mplmer setting^. The differences are that each ion needs to 
exhibit only an internal three-level structure |0), |1) and |e), and the amplitude ratio tan(6'/2) 
as well as the phase shift (j) of two laser beams should be kept constant during each pulse 
pair. The resulting Hamiltonian acting on the computational subspace {|00), |01), |10), |11)} 
reads 



uj{t) \ sm -e 


10/2 


ee)(00| 


cos -e 


*'^^^|ee)(ll| -I- h.c.^ . 


(38) 
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The TT pulse criterion oj{t)dt = vr results in the geometric two-qubit gate 

U[C] = cos0|OO)(OO|+sin0e-*‘^|OO)(ll| 

+ sin11)(001 - cos0|ll)(ll| 

+ |01)(01| + |10)(10|. (39) 

The path C, being characterized by the unit vector n = (sin 0 cos 0, sin 6^ sin 0, cos 0) in 
is traversed in the three dimensional subspace spanned by {|00), |11), |ee)} of the internal 
degrees of freedom of the ions. Note that U[C] is entangling as it cannot be written as a 
product of unitary operators acting locally on each ion qubit. It completes the universal 
set of non-adiabatic non-Abelian geometric gates that can be used to perform fast quantum 
computation by purely geometric means. 

Experiments 

The scheme for non-adiabatic GQC proposed in^ has been realized in several recent exper¬ 
iments®®. Non-commuting operations in a superconducting transmon one-qubit system 
have been demonstrated®. A universal set of one- and two-qubit non-Abelian geometric 
gates has been implemented® in a liquid state NMR quantum information processor using 
a three-qubit variant of® akin to®. 

To achieve GQC at room temperature in a naturally scalable system, spin qubits asso¬ 
ciated with a nitrogen-vacancy color centre in diamond have been used®®. While® was 
limited to one-qubit operations, the full universal set of quantum gates, including a universal 
CNOT gate, which is CNOT|a;) ® \y) = |a;) ® \xQ)y) (© is addition mod 2 and x, y take values 
0 or 1), has been implemented®. This CNOT gate applied to an initial product state has 
been shown to yield a concurrence® of 0.85, which unambiguously confirmed its entangling 
nature. 

A recent experiment® has demonstrated GQC using adiabatic evolution in a tripod 
configuration, where three ground state levels couple to an excited stated by use of three 
resonant laser fields. This experiment realized a universal set of one-qubit gates based on 
the adiabatic scheme proposed in® resulting from adiabatic transport of this dark energy 
eigensubspace for a single trapped ion. 
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Robustness of GQC 


A key motivation for GQC is the potential robustness of GPs to errors, such as decoherence 
and parameter noise^. To examine the validity of this conjecture, several studies have been 
carried out, testing the behavior of different forms of GQC to different kinds of errors. We 
summarize a few of these as follows: 


• To address the question whether geometric gates are more robust than dynamical 
gates, Zhu and Zanardi^ have designed a scheme where the two types of gates can 
be continuously changed into each other. The basic idea is to implement a one-qubit 
phase shift gates U{'j) = (-1-| -1- )(—| by exposing the qubit to a rotating 

magnetic held. By suitably changing the frequency and opening angle of the magnetic 
held, 7 can be varied from a purely dynamical phase 5 to a purely geometric phase 
<h[C]. The optimal hdelity under parameter noise is obtained for a purely geometric 
gate. Although this result provides evidence for the advantage of GQC with regard 
to resilience to dephasing errors in this setting, there are other settings where the 
geometric approach seems to have no particular advantage compared to a dynamical 
approach when considering certain decoherence^ and parameter noise^ models. 


The conditional two-qubit geometric phase gate in Eq. (29) exhibits optimal parameter 
values, where Ay is resilient to errors in the amplitude of the RF held® Similar 
optimal working points have been found at certain run-times in hnite time realization 
of non-Abelian GQC under inhuence of an oscillator bath®® ^nd parameter noise®. 


• The robustness of adiabatic and non-adiabatic realizations of non-Abelian GQC to 
parameter errors, decay, and dephasing have been examined®. The adiabatic gates 
are robust to decay and mean detuning error in the large run-time limit, but they 
are highly sensitive to dephasing and relative detuning error in this limit. The non- 
adiabatic gates, on the other hand, become resilient to all these imperfections by 
employing pulses that are sufficiently short. However, there is a limit for how short 
the pulses can be before RWA breaks down. In fact, the experiment in® did use pulse 
lengths close to this limit. Thus, one can predict that further speed-up in this setup 
would lead to highly unstable gate^®. 


21 





• Corrections to the GP of a spin-^ in the adiabatic T —)■ oo limit have been fonnd to 
decay as 1/T, while the corrections to the dynamical phase grows as T, by assnming 
a physically reasonable parameter noise model^. This effect has been stndied ex¬ 
perimentally for microwave-driven snpercondncting qnbits®®' and trapped polarized 
nltracold nentrond^. 

• Robnstness can be improved by combining GQC with other error resilient schemes. 

Adiabatic non-Abelian GQC have been combined*^ (see alsd^^^HinSl) -^ith the theory 
of decoherence free snbspaces (DFSs)®^*^ and noiseless snbsystemd^. This idea has 
been generalized to the non-adiabatic case®^. It has farther been proved that GQC 
is scalable nnder a reasonable noise model by combining it with fanlt tolerant qnantnm 
error correction^mi^. All these approaches are experimentally challenging since they 
reqnires 3- or 4-body terms in the nnderlying Hamiltonian, while only 2-body terms 
occnr natnrally in natnre. This problem may be addressed by nsing pertnrabtive 
gadgets^^ to simnlate many-body interactions. An alternative solntion is the recent 
proposal^^ to combine the Zhn-Wang approach to Abelian GQC®^ and the theory 
of DFSs, which removes the need for many-body terms bnt retains the nniversality of 
the geometric scheme. Finally, non-adiabatic GQC has been combined with dynam¬ 
ical deconpling^nS!^ j i(jga to average ont the effect of noise by hne-tnned spin 

flipping^n^. 

MIXED STATES 

The concept of pnre qnantnm states is an idealization. In experiments, mixtnres of pnre 
states arise natnrally dne to inevitable imperfections in the preparation procednre and open 
system effects dnring evolntion. Therefore, a more realistic description of experiments re¬ 
qnires the notion of mixed qnantnm states. 

The mathematical representation of mixed states is given by density operators, which are 
linear operators p satisfying p > 0 and Trp = 1. For a mixtnre of pnre states {V’j} occnring 
with relative freqnencies {pj}, one can write p = which explicitly entails the 

relation between the mixtnre and its mathematical representation p. The density 
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operator contains all empirically available information in the sense that it can be used to 
fully predict probabilities of outcomes in any conceivable experiment on a given quantum 
system. 

Now, to examine the robustness of GQC under influence of open system effects as well as 
imprecise preparation, we need to understand the meaning of GP associated with evolution 
of mixed quantum states. The basic question is, in other words, how to associate a physically 
meaningful GP to a path t G [0, r] —)■ p{t) of mixed states. For closed system evolution, such 
a path is a continuous one-parameter family of unitary transformations of the input density 
operator p(0), i.e., p{t) = U{t)p{Q)U^{t)] for an open system the evolution takes the form 
p{t) = e-^o where £ is a superoperator. For instance, in the Markovian limit case, C 

takes the Lindblad formlI^£(t)p = —{i/h)[H{t),p\ -l- {j^kpL\ ~ (1/2){P5-hl-hfc}j, where 
the Lindblad operators are arbitrary linear operators and [•, ■] ({■, ■}) is the commutator 
(anti-commutator). 

There are two main forms of GP for mixed states: 

(i) GP based on interferometry®. Here, the idea is to start from a Mach-Zehnder setup 
shown in Figure where an incoming beam carrying an internal state p (spin, say) 
is split by a 50-50 beam-splitter, the beam-pair brought together by two mirrors to 
hnally interfer at a second 50-50 beam-splitter. The internal state p is assumed to 
be unaffected by the beam-splitters and mirrors. By exposing one of the beams to a 
unitary U, acting on the internal state p, and the other beam to a U(l) phase shift e*^, 
we obtain the output intensity 

/ oc 1 -I- |Tr(pf/)| cos[x — argTr(p[/)] (40) 

in one of the output beams. Thus, the interference oscillations produced by varying 
y; is characterized by the phase shift argTr(p[/) and visibility |Tr(pl7)|. The shift 
argTr(pt/) is the Pancharatnam relative phase®22l acquired by the internal state p 
exposed to a unitary U. 

To see how this setting can be used to dehne a geometric phase associated with the 
unitary path C : t G [0,r] —)• U(t)pW(t), we use the spectral decomposition of the 
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Figure 4: Mach-Zehnder interferometer setup to measure the Pancharatnam relative phase 
argTr(pt/) acquired by an internal state p exposed to a unitary U in the upper beam. The 
phase is measured as a shift in the interference oscillations in one of the output beams 
obtained by varying the U(l) shift y in the lower beam. 

density operator, i.e., p = 'YliPn\n){n\ with being time-independent probabilities of 
the orthonormal eigenstates \n) of p. We obtain 

N N 

Tr[pl7(r)] = ^p„(u|[/(r)|n) = ^p„|(u|[/(r)|n)|e*-^T (41) 

n=l n=l 

fn being the global phase acquired by the eigenstate \n) under the evolution and 
N being the dimension of the system’s Hilbert space. Each such global phase con¬ 
tain a geometric part aig{n\U {T)\n) i {n\U\t)U{t)\n)dt and a dynamical part 
—i ^{n\U\t)U{t)\n)dt. By demanding parallel transport of each such eigenstate, we 
obtain the geometric phase shift^^ 

N 

4)[C] (42) 

n=l 

thus the GP factor is the weighted sum of GP factors of the eigenstates of p with 
weights Pn |{n|t/(r)|n)|. This GP is experimentally accessible since there are exactly 
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N independent phase factors that are not determined by the evolution, 

as can be seen by noting that the U{t)pU'^{t) is unchanged under the transformation 
U(t) —)■ These phase factors can be used to realize parallel 

transport of each eigenstate of the density operator. 


Provided p is non-degenerate, $[C] in Eq. (42) is a unique concept of GP. For a 
degenerate eigenvalue, however, the corresponding eigenstates are not uniquely given 
and therefore the above definition does not provide a unique GP. The generalization 
to degenerate density operators has therefore been carried out*^, based on the idea 
that the degenerate subspace(s) can be associated non-Abelian GPs in the same vein 
as in the Wilczek-Zee^ and Anandan^ approaches. 

The mixed state GP in^^ has been generalized and to arbitrary open system evolu- 
tion®] (for some applications, see^^^^^^^^, for which the probabilities of the spectral 
decomposition are generically time-dependent, i.e., p{t) = Yl,n=iPri{t)\'>Pn{t)){'<Pn{t)\ 
with = 5nm- The natural generalization of Eq. (42) is®^ 


N 


*|C] = arg^ Vp40KW(i/’n(0)IV'„(T))e‘J'« 


(43) 


n=l 


where the square root of the initial and final probablities is a consequence of Schmidt- 
form purification p{t) —)■ '\/Pn(^)|t/’n(^))|an), where {|a„)} is a fixed orthonormal 

ancilla basis. 


(ii) The Uhlmann This is a non-Abelian GP of mixed quantum states based on 

purifcation, i.e., the idea that any density operator can be written as the partial trace of 
a pure state of a larger system consisting of the considered system and an ancilla. This 
GP is defined via a parallelity condition that singles out a preferred purification of a 
density operator pk+i given a purification of another density operator pk. Explicitly, a 
purification of a density operator p is a mapping p —)■ \'tpp) = Ylin=ii\fP\^ri))®V^\en)■, V 
being an arbitrary unitary with T transposition with respect to the fixed orthonormal 
basis |e„). Parallelity is defined as 

max I iV'pJ I = Tr J ^/^pk+i^/^ (44) 

Tfc + l ^ 
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where the right-hand side is the Uhlmann fidelity*^ that measures the similarity of 
the two states pk and Pk+i- The optimal purifications is given by the unitary 

14+1 = (45) 

where is the unitary part of Here, Vp^,^^p^ is a unique unitary 

provided the density operators are both full rank. The lower rank cases, such as the 
pure state case, must be treated separately; in these cases the relative phases become 
partial isometries. By repeating this argument along a sequence C : pi,...,px we 
obtain the Uhlmann GP 


C'lc] = 


(46) 


where the mutliplication with U/ from the right removes the arbitrary choice of ini¬ 
tial purification, and makes the resulting unitary operator U[C] gauge invariant. The 
Uhlmann scheme can be applied to any continuous rank-preserving path by taking 


appropriate limit of the right-hand side of Eq. (46). 


One may notice that the above scheme is of rather mathematical nature and its physical 
meaning might seem a bit unclear. Nevertheless, it has been pointed out that it can 
in principle be realized interferometrically, provided the ancilla system can be fully 
controlle d^^^ * ^^^ ! Recent work have further shown the usefulness of the Uhlmann GP 
to analyze temperature effects in topological states of matter^HHlsoj 


We have seen that there are two different routes to define a GP for mixed quantum 
states. A natural question then arise: are the resulting phases in some way related? Various 
aspects of this issue have been addressed. First, it has been shown that the Uhlmann and 
interferometer GPs are indeed distinct concepts that only fully coincide in the special case of 
pure stated^^SHSH j^gy difference lies in that while the interferometer based GP can be 
measured in a single-particle interference the Uhlmann GP is a property of a larger system 
in the sense that it can only be measured in two-particle interferencd^^SIISz! Further, the 
question as to whether the two mixed state GPs arise out of a single more fundamental 
notion of GP for mixed states has been addressed^^^^®^! 
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Physical example: Unitary evolution of a qubit 

We illustrate the two GP concepts for mixed states in the case of unitary evolution of a 
single qubit. Any qubit density operators can be written as 

P=^(i + r-o-), (47) 

where 1 is the identity operator and the Bloch vector r = {x,y,z) satishes |r| < 1, with 
equality if and only if p is a pure state. Thus, the pure states reside on the two-dimensional 
Bloch sphere and the non-pure states are in the interior of this sphere. In particular, the 
origin r = 0 corresponds to a equally weigthed mixture of any orthogonal states (random 
mixture). 

Let us hrst consider the interferometric approach. Let the eigenstates of p be |±) with 
eigenvalues |(l±r), r = |r|. Assume that the qubit undergoes cyclic unitary evolution such 
that the Bloch vector r traces out a path C that encloses a solid angle G. It follows that 
|(±|G(r)|±)| = 1 and 

|Tr(pf/(r))|e*^[^l = ^ (48) 

which yields 

<h[C] = — arctan 

|Tr(pt/(r))| = ycos^ ^ r^sin^ (49) 

where the hrst expression gives the mixed state GP and the second is the visibility of the 
interference fringes that would be detected in an interferometer setting. Both these quantities 
are determined by the enclosed solid angle G on the Bloch sphere and the degree of mixing r; 
they reduce to the expected expressions <I)[C] = —f2/2 and |Tr(p[/(r))| = 1 in the pure state 
limit r —)• 1. Note that the density operator is degenerate in the limit of random mixtures 
r —)■ 0, which implies that there is no unique eigenbasis of p and the mixed state GP becomes 
undehned. Subtle interference effects close to this singularity have been examine d^^^ * ^^^ ! 

Let us now turn to the Uhlmann GP. This GP is a non-Abelian quantity, which mean 
that it involves a time ordered product or integral. To deal with this complication, one may 
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restrict to specific paths in state space. Here, we consider nnitary evolution of the qubit 
around a great circle Q inside the Bloch sphere chosen to start and end at the positive z axis 
and restricted to the xz plane. Such paths play a natural role to study temperature-driven 
phase transitions of the Uhlmann GP of fermion systems^^^MISQ] 

U[g] = (50) 

In the limit of random mixtures r —)■ 0, the Uhlmann GP is still well-defined but trivial, i.e., 
U[Q] = 1. The pure state limit r —>■ 1, on the other hand, is more subtle. In this limit, p has 
rank 1 and U[Q] should be interpreted as a U(l) phase factor times a projector. Explicitly, 
one Ends 


u\e] = -i+)(+i, (51) 

where the minus sign coincides with the expected pure state GP of tt associated with a great 
circle on the Bloch sphere. 


Measurement of mixed state GPs 


Direct experimental tests of the mixed state GP^ by using interferometry techniques have 
been carried out in NMR systems^^^^^®^. These experiments utilize two nuclear spin-| with 
one spin playing the role of the interferometer arms while implementing the unitary geo¬ 
metric transformation on the second spin, conditionally on the state of the hrst spin. The 
dependence of <h[C] in Eq. (49) on the degree of mixing r has been verihed^^^^*^, as well 
as the r and D dependence of the visibility*^. A related NMR experiment*^ has examined 
the mixed state GP proposed in®* for non-unitary evolution. 

Experimental test of*^ has been performed by using polarization mixed states of pho¬ 
tons in a Mach-Zehnder interferometer setup*®. These polarization states were prepared 
in two different ways, either by using birefringent decoherers that couple the single pho¬ 
ton’s polarization to its arrival time relative to the trigger or by preparing non-maximally 
polarization-entangled photon pairs and measuring the phase of one of the photons while 
tracing over the other. These two cases correspond to two types of mixed states, proper and 
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improper ones*^, where the latter concept refers to that the state of the full system is in 
fact pure. 

In addition to interferometry, quantum mechanical phase shifts can as well be performed 
using single-beam polarimetric techniques. The basic setting to perform such experiments 
has been developed® and later generalized to mixed states^^®. The latter has been realized 
in a single-beam experiment on partially polarized neutrons!^ and utilized to test the non¬ 
additive nature of phase shifts (geometric or dynamical) for mixed quantum states^^. 

Due to the high-level of control required, the Uhlmann GP is considerably more chal¬ 
lenging to implement experimentally. Nevertheless, it has been measured by using NMR 
technique for a qubit undergoing unitary evolution and comparing it with the interfero¬ 
metric GP for the same paths in the space of density operators*^ . The measured phases 
were explicitly shown to behave differently, providing clear experimental evidence of the 
inequivalence of the interferometer-based and Uhlmann GPs. 

ENTANGLEMENT 

Entanglement is fundamental resource in quantum information that can be used for secure 
key distribution*^, teleportation*^, and information processing®. The most basic form of 
entanglement is given by the four Bell states, which form the orthogonal basis states 

| 4 .±> = ^(| 00 )±| 11 )), 

l>t±> = ^(|01)±|10)) (52) 

of two qubits. These are maximally entangled in the sense that each of them contains 
full information about the two qubits, but completely random information for each qubit 
separately. 

The concept of entanglement becomes even richer when considering more than two qubits. 
The additional richness lies in the fact that while two qubits essentially can be entangled in 
only one way (given by the above Bell state entanglement), there are several ways to entan¬ 
gle three or more qubits. To analyze this structure, it is convenient to use the concept of 
stochastic local operations and classical communication (SLOGG), i.e., the idea that entan- 
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glement cannot be changed if only local operations (nnitary, measurements, etc) and classical 
communication are allowed. In this language, one may for instance show that three-qubit 
system in pure states can entangle in two different ways, in the sense that a state picked 
from one of these two classes cannot be taken to a state in the other class by performing 
SLOCC^^. The two classes, called W and GHZ, have their basic forms 

\W) = a|001)+6|010)+c|100), 

|GHZ) = a'lOOO)+ 6'|111), (53) 

where the hrst is characterized by two-qubit entanglement only, while the second is genuinely 
three-qubit entangled but with no two-qubit ent anglement*^. 

Entanglement can be quantihed by using various kinds of measures, such as entangle¬ 
ment of format ion^^, geometric entanglement^^^^nSSl, ^nd relative entropy of ent anglement*^. 
Besides quantihcation of its amount, entanglement can be characterized qualitatively by 
using suitable polynomials of the expansion coefficients with respect to a product basis 
spanning the composite system’s Hilbert space. For two qubits in the state 'Yliki^ki\kl)., 
det[afci] = aoo«ii — «oi«io is a polynomial invariant under SLOGG and therefore a property 
of entanglement. Furthermore, this degree-2 polynomial determines the entanglement of 
formation via the concurrence measure^ C = 2 |det[Q!fci]|. Similarity, the hyperdeterminant 
of the expansion coefficients akim of a pure three-qubit state is a SLOGG invariant degree-4 
polynomial that determines the 3-tangle, one of the main measures of three-qubit entan- 
glementi^. The characterization in terms of invariant polynomials becomes increasingly 
complicated when the number of subsystems increases^^^^^^. 

In the following, we focus on a type of topological phases for entangled systems undergoing 
local special unitary (SU) evolution; a subclass of SLOGG operations. These topological 
phases were hrst discovered for two-qubit systems^^^^^^^^"^^, and later generalized to qudit (d- 
level) pairs®^^®^^ and to multi-qubit system^^. They may in some cases coincide with the 
corresponding GPs. We delineate the relation between the topological phases and local SU 
invariant polynomials and describe experimental work to study these phases. The key merit 
of the entanglement-induced topological phases is that they constitute a novel perspective 
that may provide further insights into the nature of quantum entanglement. 
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A non-exhaustive list of other studies of GP effects in entangled systems is as follows. GPs 
of bipartite qubit systems in pure states have been examined*^^^*^^ in terms of the concept of 
‘Schmidt sphere’^^, built on the analogy of the Bloch sphere of a single qubit, and by use of 
properties of braiding transformations to represent entangled two-qubit stated^. Recently, 
the Schmidt sphere concept has been examined in a quantum optical experiment*^. The 
relation between the GPs of an entangled system and its subsystems has been examined* ^^^ * ^^^ *. 
GP effects have been shown to play a role for the understanding of entanglement in multi¬ 
particle systems*^^^^'*^. GPs for sequences of relative states, obtained by projecting on one 
subsystem of a bipartite composite system, have been studied for pure*^ and mixed*^ 
states. GPs associated with spin-orbit entanglement in neutron interferometry have been 
predicted*^ and experimentally observed*^. A relation between the GP of a two-qubit state 
and concurrence has been delineated*^. Pure two-qubit states can be represented as single¬ 
qubit states with quaternionic expansion coefficient^^. The non-Abelian and entanglement- 
dependent GPs associated with this quaternionic representation have been examined* ^^^ * ^^^ *. 


Topological phases of entangled systems: bipartite case 

Let us hrst consider the case of maximally entangled states (MES) of pairs of qubits. Any 
such state can be written on the form 


a,b) = i(^a|00) + ()|01> - r|10) + a*|ll)j, 


(54) 


where a and b are arbitrary complex numbers with |ap -|- |6p = 1. The four Bell states in 


Eq. (52) correspond to the choices (a, 6) = (1,0), (i,0), (0,i), and (0,1), respectively. These 
states can be transformed into each other by applying local SU(2) transformations U and 
V. The states of the single qubit subsystems A and B are given by the reduced density 
operators pA and ps obtained by tracing over qubit B and A, respectively. For the MES 
|a, b), the reduced density operators take the form and \1b, irrespective of the values of 
a and b. Thus, the subsystem states correspond to random mixtures for all MES. 

We are interested in continuous paths C : t G [0,r] — )■ {a(t),b{t)) corresponding to 
the pure state local unitary evolution |a(0),6(0)) — )■ \a(t),b(t)) = t/(t) ® R(t)|a(0), 6(0)), 
such that the state traces out a loop in state space of the qubit-pair, i.e., |a(r),6(r)) = 
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e*-^|a(0), 6(0)). We can now observe two important features of this evolution; 


In order to preserve the MES form in Eq. (54), / can only take one of two values, 
namely, 0 or tt. 


The accumulation of local phase changes 
i [ {a{t),b{t)\ — \a{t),b{t))dt = i [ Tr () 


dt 


dt + i Tr -1 






LS' 




dt (55) 


vanishes, since the inhnitesimal generators UW and EW of U, V gSU( 2) are traceless. 


The hrst feature demonstrates the topological nature of the acquired phase /: it is 
impossible to transform between the two allowed values 0 and tt by continuously deforming 
the path C. In fact, the two phases can be associated with the doubly-connectedness of 
SO(3). This may be seen by noting that the state space of MES is given by 

MES ~ {(a, b) G such that |ap + |6p = 1 and (a, b) ~ (—a, —b)}. (56) 

This shows that there is a one-to-one correspondence between MES and the three-dimensional 
sphere (|ap -|- |6p = 1) with antipodal points identihed ((a, 6) ~ (—a, —6)), which is iso¬ 
morphic to the three-dimensional rotation group SO(3), i.e., S^/Z 2 — SO(3). Now, SO(3) 
is known to be doubly-connected, which means that there are two types of topologically 
inequivalent loops, one which is trivial, i.e., can be continuously deformed to a point, and 
one which is trivial only when traversed twice. These two types of loops are illustrated in 
Figure and correspond to the two topological phase factors -|-1 and —1, respectively, first 
discovered by Milman and Mosseri^. 

The second feature demonstrates that / is the GP associated with the loop since the 
accumulation of local phase changes along the path vanishes. However, the local phase 
changes vanish only for maximally entangled states; for more general two-qubit states It/j) = 
a|00) -|- 6|01) -)- c|10) -|- d|ll), the global phase factors contain also dynamical contributions, 

i.e., 

i j {'ip\'ip)dt ~ ^ y {^PaU^U^ dt + i J Tr dt ^ 0, (57) 

since pa and Pb 7 ^ if IV’) is non-maximally entangled, i.e., if 2 \ad — bc\ < 1. 
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Figure 5: Illustrating the two classes of topologically inequivalent loops in S'^/Z 2 — SO(3). 
The loop C (solid line) can be continuously deformed to a point and correspond to the trivial 
phase factor +1. On the other hand, C (dashed line) can be continuously deformed to a 
point only when traversed twice and correspond to the non-trivial phase factor —1. Note 
that C is a loop since antipodal points on the sphere are identihed, i.e., P and Q are actually 
the same point. 

Note that the allowed cyclic phases / are still 0 or tt for any non-vanishing degree of 
entanglement. The only exception is when l"^) is a product state, i.e., when ad — be = 0, in 
case of which a continuum of cyclic phases / may occur under local SU(2) evolution. Thus, 
the topological nature of the allowed global phases / is intimately related to entanglement, 
in the sense that / is no longer a discrete quantity exactly when 1-0) ceases to be entangled. 
Note that the degree-2 polynomial ad — be is the determinant of the expansion coefficient, 
which, as noted above, is known to be a SLOOC and thereby local SU invariant. 

The topological two-qubit phases have been generalized by Oxman and Khoury to pairs 
of d-level systems (qudits)®! and pairs of systems with different dimension*^. For equal 
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dimension, they found the allowed phase factors to be the dth roots of unity 




(58) 


and being related to topological properties of SU{d). The above result for qubit-pairs is 
recovered as the special case where d = 2. 


The discrete set of phases in Eq. (58) can be understood as follows. Let the initial 


state take the general form |'0(O)) = Ylti=i ^ki\kl). Define the d x d matrix ct with compo¬ 
nents A local SU(d) evolution |'0(O)) —t U{t) ^ l/(f)|'0(O)) can be translated into the 
matrix evolution cx i—>■ u(t)a\"^(t), where v(f) and u(t) are SU(d) matrices with elements 
Uki{t) = {k\U{t)\l) and Vki{t) = {k\V{t)\l). Now, a cyclic evolution |'0('r)) = yields 

u(r)Q:v"'"(r) = By taking the determinant of this expression, we obtain 


det[u(r)] det[Q:] det[v"'"(r)] = det[Q:] = det[Q:], 


(59) 


where we have used that det[u(r)] = det[v'^(r)] = 1 for SU{d) matrices. If det[Q:] 7 ^ 0, we 


thus obtain that = 1, which implies Eq. (58). On the other hand, if det[Q:] = 0, which 


happens if and only if |' 0 (O)) is a product state, then / can take any value and is thereby 
no longer topological. Note that this is consistent with the result we found in the two-qubit 
case discussed above. 


Topological phases of entangled systems: multi-qubit case 

The topological phase structure of multi-qubit systems undergoing local SU(2) evolution 
has been examined^ZI. These phases show a considerably richer structure than those in the 
bipartite case, due to the richer entanglement structure of systems consisting of more than 
two particles. 

As an illustration of this additional richness, let us consider the simplest case of three 
qubits in some detail. Here, the allowed topological phases are 0, |, 7 r, To under¬ 
stand the reason why just these phases, let us consider an arbitrary three-qubit pure state 
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J2kim'^kim\klm). The hyperdeterminant of the expansion coefficients akim takes the form 

Det[afc;m] = Q^ooo'^iii + '^ooi'^iio + '^oio'^ioi + '^lOO'^Oll 

~2 (aoooctooictiioctiii + ctoooctoioctioictiii + ctoooctiooctoiictiii 
+aooiQ!oioctioiQ!iio + ctooictiooctoiictiio + ctoioctiooc^oiictioi) 

+4 (aooo«oii«ioi'^iio + «ooiOoio«ioo«iii), (60) 

which apparently is a degree-4 polynomial. Just as det[afc;] in the bipartite case, Det)^^;^] 
is a SLOCC invariant; in fact, it determines the 3-tangle = 4|Det[aMm]|- The 3-tangle 
measures the amount of genuine three-qubit entanglement in the state*^. 

Now, given a cyclic local unitary evolution where all akim —>■ e'^^akim, we obtain Det[afcZm] —t 
Det[afc;m] = Det[afc;m] by using the local SU invariance of Det)^^;^]- Thus, provided 

Det[aMm] 7^ 0, we find = 1, which implies 

/ = g(27r/4), g = 0,...,3. (61) 

As in the bipartite case, we see that / can take any value and whereby losing its topological 
nature, in cases where 'Dei[akim\ = 0, which is exactly when three-qubit entanglement is lost 
(i.e., Ts = 0). In this sense, the topological phase / is a genuine property of three-qubit 
entanglement. 

By looking at the space of all pure inseparable three-qubit states, we note that the In¬ 
states is a zero measure set of this space with vanishing 3-tangle. The IT-set constitutes a 
singularity in the space of entangled three qubits, just as the singularity line defined by the 
magnetic solenoid that gives rise to the Aharonov-Bohm topological phase*^. This further 
illustrates the toplogical nature of the allowed / for non-vanishing 3-tangle. 

The number of allowed topological phases increases with the number of qubits. They 
have been found for up to seven qubits by using a search algorithm based on a combinatorial 
formulation of the problem^. For eight or more qubits, only partial results have been estab¬ 
lished due to rapid increase in the computational resources needed to execute the algorithm 
when the number of qubits increases. 

Finally, we note that the multi-qubit phases may coincide with the corresponding GPs. 
As in the bipartite case, this happens precisely when all the marginal one-qubit states are 
random mixtures. 
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Experiments on topological phases of entangled systems 

Souza et (see alscA^ have used an optical setup involving non-separable polarization 

and orbital degrees of freedom to simulate quantum entanglement in a laser beam. The 
laser beam states relevant for this experiment can be described by multiplying orthogonal 
polarization unit vectors and sy (horizontal and vertical polarization, respectively) with 
a pair of first order Laguerre-Gaussian profiles t/’±(r) to define the laser beam amplitude 

E(r) = ^(^a^/J+{r)eH + I3ip+{r)ev - l3*ilj-{r)eH + a*^/J-{r)ev^ (62) 

representing an arbitrary maximally non-separable (MNS) mode, being a laser beam analog 
of MES for two-qubit systems. Here, a, (3 are complex numbers such that |ap -f |/3p = 1. 
All these states are locally equivalent in the sense that any MNS mode can be reached 
by manipulating only the polarization, say, of the beam. This fact was utilized in*^ to 
implement loops in the space of MNS modes by letting the laser beam pass a sequence of wave 
plates. In this way, the topological phases associated with these polarization transformations 
were realized and observed. 

An analogous experiment to observe the topological two-qubit phase for MES has been 
performed by using NMR techniqud^. The setup consists of two nuclear spin-qubits pre¬ 
pared in a MES and a third ancillary qubit playing the role of the two arms of a Mach-Zehnder 
type interferometer. The input state :^(|0) -|- |1)) ® |MES) is acted on by a conditional op¬ 
eration that takes only the MES copy connected to the |1) state around a loop. Thus, the 
three qubit state undergoes the transformation 

^(|0) + |1)) 0 IMES) ^ ^(|0> ± |1)) 0 |MES), (63) 

where the relative sign is precisely the topological phase factor ±1 associated with the loop. 
In this way, the topological phases were read out by measuring the ancilla qubit state*^. 

The qudit and multi-qubit phases have not been verified experimentally yet. However, 
feasible experimental proposals using polarization and orbital degrees of freedom of photon 
pairs have been put forward^^^^^^ and awaits to be performed in the future. 
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CONCLUSIONS 


We have reviewed the role of geometric phase (GP) ideas in quantum information science. 
We have described how GPs can be used to perform robust quantum information processing, 
leading to the idea of holonomic or geometric quantum computation, which is a tool to realize 
quantum gates that are robust to certain errors. GP concepts for mixed quantum states have 
been described, focusing on Uhlmann’s seminal work from the mid 1980s and the more recent 
development of mixed state GPs in interferometry. Finally, we have decribed how GPs can 
be used to analyze quantum entanglement, with particular focus on the discovery of a new 
type of topological phases for entangled systems. 

The importance of the reviewed body of work is two-fold. It has led to new ways to 
perform quantum information processing that may be useful to hght errors in quantum gate 
operations. This may help to reach below the error threshold, below which quantum error 
correction codes can be performed. In other words, GPs may help to realize large-scale 
quantum computation. Gonversely, ideas that have been developed or rehned within the 
quantum information community, such as the theory of mixed quantum states, open system 
effects, and quantum entanglement, have been applied to the theory of GP. In this way, 
further insights into the physical, mathematical, and conceptual nature of the GP have been 
obtained. 

There are a number of pertinent issues to examine in the future in this research area. 
First, explicit physical implementations of schemes that combine GQG and symmetry-aided 
error protection techniques such as decoherence free subspaces and subsystems, need to be 
developed and experimentally implemented. The key question here is whether there exist 
such schemes that at the same time avoid many-body interaction terms in the Hamiltonian 
and thereby would be much simpler to realize in the laboratory than the existing schemes. 
Secondly, the conjectured robustness of GQG compared to more conventional dynamical 
quantum gates is still an open problem that needs to be examined further. In particular, 
the robustness features of geometric gates in the presence of open system effects need to 
be addressed in realistic models that include signihcant non-Markovian effects. Thirdly, the 
relation between the interferometric-based and Uhlmann GPs is still not fully understood. A 
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key goal here is to find a common conceptual basis for these two phases, such as for instance 
a new mixed state GP concept that covers the two. Finally, a major challenge concerning 
GP effects in entangled systems is to clarify the relation between the entanglement-induced 
topological phases and SLOGG invariant polynomials of multi-qubit systems and to imple¬ 
ment experimentally these phases in the qudit and multi-qubit cases. Addressing the above 
listed open problems would lead to further insights into the nature of the GP and to new 
applications of GP ideas in quantum information science. 
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